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SUMMARY 


This report compares the simulation results of a real time control algorithm for onboard 
computation of time-optimal intercept trajectories for the F-8 and F-15 aircraft. Due to the inherent 
aerodynamic and propulsion differences in the aircraft, there are major differences in their optimal 
trajectories. The significant difference between the two aircraft are their flight envelopes. The F- 
8's optimal cruise velocity is thrust limited, while the F-15's optimal cruise velocity is at the 
intersection of the Mach and dynamic pressure constraint boundaries. This inherent difference 
necessitated the development of a proportional thrust controller for use as the F-15 approaches it's 
optimal cruise energy. Another interesting phenomena is that the optimal climb trajectory for the F- 
15 is along its dynamic pressure boundary. This necessitated the use of a sub-optimal proportional 
vertical lift controller to track the constraint boundary. 

This report documents the application of singular perturbation theory to the trajectory 
optimization problem, along with a summary of the control algorithms. Numerical results for the 
two aircraft are compared to illustrate the performance of the minimum time algorithm, and to 
compute the resulting flight paths. A major recommendation is that future research be directed at 
the application of singular perturbation methods to problems in flight mechanics where state 
constraints, such as a maximum dynamic pressure limit, play an important role in the analysis. 

This report documents a portion of the total research effort supported under this grant. The 
research related to time optimal aircraft pursuit evasion can be found in [4]. 
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SECTION 1 
INTRODUCTION 


There has been active research into optimizing flight path trajectories using multiple time 
scale separation techniques. In [1] Calise and Meorder showed that by using singular perturbation 
theory to separate the dynamics into fast and slow modes and then applying the optimality 
conditions from calculus of variations, one could obtain a closed form solution to the min-time 
intercept problem. What soon followed was a 3-D real-time piloted simulation and which used 
aerodynamic and propulsion data from the F-8 [2]. The algorithm was eventually flight tested on 
the NASA F-8 test aircraft at NASA Ames/Dryden Flight Research Facility [3]. The objective of 
this research was to modify the existing algorithm for use on an F-15 aircraft. 

This report documents the results of modifying the min-time intercept algorithm from the F- 
8 to F-15, which represents a portion of the activity supported under this research grant. A portion 
of this work was conducted at NASA Ames/Dryden Flight Research Facility during the summer 
1989. This report will contain four sections. The problem formulation is given in Section 2. A 
summary of the control algorithm is given in Section 3, which highlights the application of singular 
perturbation theory along with the optimality conditions derived from calculus of variations. This 
development could have been explained in full detail, however, in the interest of brevity only the 
first two boundary layer approximations will be discussed. This also coincides with the fact that a 
sub-optimal proportion lift controller was used on the F-15 due to the fact that the optimal 
trajectory rides the dynamic pressure constraint. Section 4 presents numerical results comparing 
the F-8 and F-15. Section 5 gives the conclusions for this work and identifies future work which 
needs to be accomplished prior to a flight test. 
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SECTION 2 

PROBLEM FORMULATION 


The point mass equations of motion are referenced to a horizontal, target centered, inertial 
coordinate frame illustrated in Fig. 1 


x = Vcos y cos P 

(1) 

y = Vcos y sin P - Vycos yp 

(2) 

eE = (T-D)V/W 

(3) 

e 2 P = Lsin |i / mVcos y 

(4) 

h = Vsiny 

(5) 

y = (Lcos p - Wcos y) / mV 

(6) 


These equations are valid for constant weight, thrust aligned with the flight path, and flat earth 
approximations. E = h + V 2 / 2g is the total aircraft energy per unit weight, P is the heading angle, 
h the altitude, y the flight path angle, and (I the bank angle. Drag is assumed to have conventional 
parabolic form 

D = qSCD 0 + KL 2 /qS, q = pV 2 /2 (7) 

where q is the dynamic pressure, p the air density and 

K = ti/CL a (8) 


Lift is defined by 


L = qSCL a a 


(9) 


where a is the angle of attack. The control variables are L, |i, and thrust T. The objective is to 
find the controls L, ji, T that minimize the time to intercept a constant velocity target 



( 10 ) 


The minimization is subject to the following state and control variable constraints: 
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L < W Gmax 
L < qS CL a CXmax 
Tmin(h,V) < Tmax(h,V) 
q ^ qmax, V < V max 


(ID 

( 12 ) 

(13) 

(14) 


where Gmax is the maximum load factor, (Xmax the maximum angle of attack, Tmin and Tmax the 
thrust level limits that are functions of aircraft altitude and velocity. The boundary conditions are 
such that the initial aircraft state is fully specified and require 


x(tf) = y(tf) = 0, h(tf) = hj(tf) 

for intercept, where hx(tf) is taken as the projected target motion in altitude 

hj(tf) = hx(0) + (VxsinYr)tf 


(15) 


(16) 


The parameter e designates multiple time scaling used to order the dynamics [1]. The approach 
here is to find a solution for e = 1 by an power series expansion around e = 0. The boundary 
layers are separated by rescaling time as = t / £*, i = respectively, then setting e = 0 in the 

resulting equations. A justification for this specific ordering of the dynamic equations is given in 

[ 1 ]. 
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SECTION 3 

SUMMARY OF CONTROL ALGORITHM 


3.1 Outer Solution 

In the outer solution, the controlled aircraft is assumed to be traveling on a fixed course at a 
constant speed, as can be seen by letting e-»0 in Eqs. (3-6). This addresses only the x and y 
dynamics, and the states (3, h, and E take on the role of control like variables. In order to satisfy 
the intercept requirements, we have the following constraint 

Vsin (p-X) = V T cos rr cos X (17) 

This means that there is no relative motion allowed perpendicular to the horizontal line-of-sight 
axis. The optimal controls ho and Eo are determined from minimizing the reduced Hamiltonian 

Ho (E,hi) = X x o Vcos (3 + ^.yo (V sin p - Vt cos yr) +1=0 (18) 

and it is shown in [1] that this reduces to 

ho, Eo = arg max (V) (19) 

h,E 

where the maximization takes place subject to the constraints in Eqs. (11-14) and subject to the 
following conditions that result from setting e=0 in (3-6) 

T = D 0 , Mo = 0, yo = 0, Lo = W (20) 

The term Do in (20) is drag for L = W 

D 0 = qoSCD 0 + KW 2 /q 0 S (21) 

where 

qo = p(h 0 )V 0 2 / 2, V 0 = V2g(E 0 -h 0 ) (22) 

The subscript 0 denotes the zeroth or outer solution. The maximization in (19) is equivalent to 
finding the maximum velocity cruise point. 

The cruise point solution from Eq. (19) for the F-8 and F-15 aircraft are displayed in Figs. 
2 and 3, superimposed upon their flight envelopes. Note that since the F-8 is thrust limited, its 
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cruise point lies on the T m ax=D contour. The F-15 is q and Mach limited, and its cruise point lies 
at the intersection of these constraint boundaries. Thus T does not equal T max in the outer solution 
for the F-15. 

The optimal cruise heading Po is computed using (17) 

Po = sin' 1 (Vj cos Yr cos Wo) + X (23) 

The costates Xxo and Xyo, associated with the horizontal position dynamics in the outer control 
solution are needed in subsequent boundary layer solutions. These take the form 

Xxo = ' cos po / (Vo - Vt cos Yr cos po) (24) 

A,yo = -sin Po / (Vo - Vt cos Yr cos Po) (25) 

These are determined from the optimality conditions 5Ho / dp = 0 and the condition Ho = 0. It 
should be noted that the cruise solution for ho and Eo is independent of target motion and intercept 
geometry. This allows these quantities to be calculated off line and stored. The only outer solution 
calculations performed on-line are Eqs. (23), (24), and (25). 

3.2 First Boundary Layer Solution 

The first boundary layer addresses only the energy dynamics. The constraints 

pi = 0, Yl=0, Li=W (26) 

arise from Eqs. (4-6) when the time transformation x = t / e is introduced and we let e->0. The 

control like variables are T, h and p. The optimal P is the same as in the outer solution. Since T 

appears linearly in the Hamiltonian, 

Ti = T max (h,V), when XEi < 0 (27) 

Ti = T min (h,V), when >.Ei > 0 (28) 

This corresponds roughly to an energy climb and energy descent, respectively. Optimization with 
respect to h yields for climb 


5 



( 29 ) 


, c f (Tmax- Do)V 1 

hi - arg max | Y y G J E = Ecurrent Tmax > D 0 

and descent 

,d - f (Tmin- Do)V 1 

hr - arg min y J E = Ecurrent Tmax < D 0 

h V-Vo 

The climb path to cruise for the F-8 and F-15 are superimposed on the aircraft flight 
envelopes in Figs. 2 and 3. The optimal descent path for both aircraft is along the q ma x boundary. 
The expression for the first boundary layer costate is 

XEi = -WHo(E,hi)/Vi(Ti-D 0 ) (31) 

where Flo is the outer solution Hamiltonian evaluated at the first boundary layer conditions. Since 
the solutions for hi c and hi^ are independent of target motion, they can be precomputed and stored 
as a function of E. Only the energy costate variable in (31) is computed on-line. 

3.3 Second Boundary Layer Solution 

The second Boundary layer solution determines the optimal heading angle dynamics. 
Introducing the time stretching transformation x = t / and letting £— >0 while holding the slow 

dynamics x, y and E fixed, yields the constraints 

Y 2 = 0, L2 = L2 22 + W2 (32) 

where L is the total lift and L 22 = Lsin p., the horizontal component of lift. The control variables 
are T, h and L 22 . Assuming that all turning takes place near the initial time where XEi < 0, the 
optimal thrust is 

T 2 = T max (h 2 ,V 2 ) (33) 

where h 2 is the optimal commanded altitude determined by 

h 2 = arg min {-p / KV Hj (E, h, P)}l E=Ecurrent, (3=Pcurrent 
h 

In Eq. (34) H^E, h, (3) is the Hamiltonian in the first boundary layer evaluated at the current 
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values of E, h and (5. It is expressed as 

Hi(E, h, P) = taco Vcos P +^.yo (Vsin P - Vj cos yr) + X,Ei [(T - Do)V / W] + 1 = 0 (35) 

The solution for the horizontal lift component is 

L 22 = V-qSWHi(E,H,|3)/(VKXf)*sign(fVp) (36) 


3.4 Proportional Vertical Lift Control 


An option was included in the algorithm for stopping the singular perturbation analysis 
after the second boundary layer, and employing a sub-optimal vertical lift solution. It was 
necessary in this study to use this option due to the fact that zero order singular perturbation 
analysis results in a steady state error when following ramp like altitude commands. Although this 
steady error is not a critical factor for the F-8 study in [1], it is essential to accurately track the 
altitude command for the F-15 since the optimal climb path lies essentially along the dynamic 
pressure constraint (Fig. 3). The derivation of the control logic proceeds as follows. From Eq. 
(5) 


h = Vsin y 


(37) 


We would like the altitude rate to track the altitude error according to proportional feedback control 

x h h = h 2 -h. (38) 


where Tj, is the time constant associated with the decay in altitude error. Since h 2 is not a constant 
command signal, but more like a ramp command, a term must be included to account for the steady 
state error that would otherwise result. The rate of change in I 12 can be estimated using 


h 2 = hi = — E = Visin Yi 
3E 


(39) 


where 3hi/9E is the slope associated with altitude with respect to E along the dynamic pressure 
boundary (Fig. 3). In (39) we have used the fact that h 2 (E) = h x (E) when the heading error is near 
zero. The desired flight path angle Yd is formed by summing y from (37) and (38) with Yi, from 
(39). Using small angle approximation for sinywe have 
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_ E 9h | (h2-h) 
Vi9E ThV 


( 40 ) 


With Yd defined we can calculate the vertical lift component. Using Eq. (6) 

Y = (Li- Wcos y) / mV = Xy (Yd - Y) (41) 

where Xy is the time constant associated with the decay in flight path angle error. As an alternative, 
Xyand Xh can be related to a desired natural frequency and damping ratio [1]. Solving for Li 
yields 

Li = mVXy (Yd - Y) + Wcos y (42) 

The relationship between the vertical and horizontal lift components along with the bank angle are 
given by 

p = arctan ^ L = V L? + L 2 


where L -2 is the horizontal component of lift. 

3.5 Proportional Cruise Thrust Control 

Since the flight envelope for the F-8 has as its V ma x cruise point at a T=D point the throttle 
setting was set at full throttle. This means that the F-8 will asymptotically reach its optimal cruise 
energy. In other words it never theoretically attains the desired optimal cruise energy. The F-15 
flight envelope has as its Vmax cruise point at the intersection of the qmax limit and M m ax limit. 
With the cruise point at this location the F-15 reaches its optimal cruise energy in finite time. By 
virtue of its high T/W ratio the F-15 can gain energy at a much higher rate compared to that of the 
F-8. Therefore in the case of the F-15 it is necessary to throttle back as the cruise energy is 
approached. This throttling was incorporated into the algorithm by use of a proportional 
controller, the derivation of which is described as follows. Since it is desired that the energy rate 
track the error in energy, which is constant during cruise, we have from Eq. (3) 

E = (T-D) V/W = (Ec-E)/xe (44) 


Solving for T gives 
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T = W(E c -E)/t e V + D 


( 45 ) 


where x E is the controller time constant set to give the desired rate of decay in the energy error. 

The procedure for throttle control during descent was the same as that used in the F-8 study 
described in [1] 
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SECTION 4 

NUMERICAL RESULTS 


Figures 4 and 5 show the ground tracks for Cases 1-3. Both simulations have the same 
initial geometry. The target flies at a constant altitude, velocity and heading but with differing 
downrange initial distances. These differing initial distances demonstrate the algorithms range- 
matching calculation of an optimal pseudo-cruise point [1]. Notice the rapid closure rate of the F- 
15 as compared to the F-8. Figures 6 and 7 show the altitude and optimal commanded altitude (h 2 ) 
time histories. The time constants Th and Ty were set to correspond to a damping ratio and natural 
frequency of 1.0 and 0.1 respectively. Case 1 is a long range intercept which is comprised of three 
phases; an initial climb along the optimal climb path followed by a cruise leg at the aircraft's 
optimal cruise energy (E*, h*), then a descent portion along the dynamic pressure boundary to 
intercept the target. Cases 2-3 being shorter intercepts never obtain their optimal cruise energies 
due to their proximity to their targets. These cases obtain pseudo-cruise energies instead [1]. In 
these cases the aircraft climbs to a pseudo-cruise energy and immediately initiates descent. Notice 
in each of the F-8 simulations the optimal commanded altitude has peaks between 150 and 175 
seconds for each case. These demonstrate the characteristic dive in the transonic region. In 
contrast, the F-15 has enough thrust to simply ride the dynamic pressure boundary. Another 
interesting feature can be seen in comparing the Case 2 time histories. The F-8 never reaches its 
optimal cruise energy. Instead it approaches a lower pseudo-cruise energy. The F-15 has 
sufficient energy rate capacity to reach its optimal cruise energy at this shorter intercept range. The 
Case 2 optimal climb paths are compared in Figs. 8 and 9. In the F-8 simulation notice the near 
constant velocity energy climb followed by the characteristic transonic dive. The F-15 climb 
profile is along the dynamic pressure boundary. There exists a slight dynamic pressure boundary 
violation in the figure, however, this can be improved by tuning the gains th, Ty and te 

Figure 10 is a comparison of optimizing the h and y dynamics [1] versus using the sub- 
optimal proportional lift controller. Notice the underdamped response of the optimized boundary 
layers results in severe dynamic pressure constraint boundary violations. This is due to separating 
the h and y dynamics into separate boundary layers which has to be done in order to obtain a closed 
form solution. The suboptimal proportional lift controller allows the designer to pick the gains so 
as to avoid a qmax violation. Figures 1 1 and 12 show the flight path angle and desired flight path 
angle time histories. In both simulations, descent initiation is clearly evident by the sudden 
decrease in the desired flight path angle. Also notice that the flight path angle never exceeds 0.2 
rad during climb, which validates the y = 0 approximation in the first and second boundary layer 
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analysis. It can also be seen that the F- 15 completes the intercept in nearly half the time it takes the 
F-8. The lift and bank angle profiles are shown in Figs. 13 and 14. Note that both aircraft bank 
180 degrees in order to initiate descent. In the F-15 simulation, the sudden decrease in lift at 1 15 
seconds is due to the approaching optimal cruise velocity, where the throttle control in (45) is 
initiated. This changes the energy rate term in (40), and consequently the vertical lift calculated in 
(42). In Figs. 15 and 16 the thrust histories for each aircraft are given. The time constant te here 
was set to 6 seconds. Descent initiation is evident by the sharp decrease in thrust at 150 and 350 
seconds for the F-15 and F-8 respectively. In the F-15 simulation, at 115 seconds the proportional 
thrust controller throttles back the engines just prior to the aircraft obtaining its pseudo-cruise 
energy level. 

The ground tracks for Case 4 are given in Figs. 17 and 18. This case is a close in intercept 
which has the target moving at constant altitude, velocity and heading 180 degrees in the opposite 
direction. This case illustrates a combined initial hard turn and climb (yo-yo maneuver) followed 
by a descent under near tail chase conditions. The pursuer's initial altitude is identical to that of the 
target and the velocity is above the comer velocity for that altitude. Note that the downrange 
intercept distance is considerably less for the F-15 aircraft. Figures 19 and 20 give the altitude and 
commanded altitude time histories. Both aircraft perform a high speed yo-yo maneuver in order to 
trade speed for increased turning performance. It is interesting to note that the time to complete the 
hard turn is nearly identical for both aircraft, 10 seconds. The reason for the large altitude 
command is due to the large initial heading errors. The jump in commanded altitude is at the 
completion of the hard turn so as to get on the optimal climb path. Between 10 and 20 seconds the 
F-8 dives to trade potential energy for kinetic energy, while the F-15 just accelerates to intercept. 
Figures 21 and 22 show altitude versus velocity plots for Case 4 where the high speed yo-yo 
maneuvers are more evident. Both simulations show the aircraft moving initially toward the comer 
velocity to trade speed for increased turning performance. Once the turn is completed both aircraft 
move to get on the optimal climb path which is followed by a descent along the dynamic pressure 
boundary. Figures 23 and 24 show the lift and bank angle profiles. Both aircraft initiate 
maximum G turns at 0-10 seconds to reduce the heading error. This is followed by a climb phase 
followed by cruise then descent at 40 and 90 seconds for the F-15 and F-8 respectively. In the F- 
15 simulation it is evident that throttling is taking place at 28 seconds due to the decrease in lift. 
Also note the large lift at interception, this indicates that the F-15 is tending to miss the target. This 
is due to the higher degree of coupling in the dynamics for the F-15, which requires further 
investigation. 
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SECTION 5 

CONCLUSIONS AND FUTURE RESEARCH 


The differences in the numerical results between the F-8 and F- 15 trajectories lie mainly in 
each aircraft's ability to gain energy. The F-8 is an aircraft which never encounters any state 
constraint boundaries during its climb phase. The F-15’s higher thrust to weight ratio allows the 
aircraft to gain energy at a much higher rate. Since the F-15's flight trajectory lies on the dynamic 
pressure boundary, the methodology for optimizing the altitude and flight path angle dynamics in 
[1] does not apply. This was circumvented by the application of the sub-optimal proportional lift 
and thrust controllers which provided an adequate sub-optimal solution. It was noticed that setting 
the time constants Xh, Xy and te to achieve the best performance for long range intercepts gave poor 
close-in intercept performance. This problem would imply that these gains are maneuver 
dependent, which it not satisfactory for real time implementation. A solution to this problem is to 
attempt to optimize the altitude and flight path angle dynamics, subject to the dynamic pressure 
constraint. 

The major recommendation for future research is that the singular perturbation 
methodology in [1] be extended to address the issue of state constrained optimization problems. 
This would avoid the gain scheduling issue described above. Another interesting point worth 
investigating is the assumption of thrust aligned along the velocity vector. Thrust is actually 
aligned with the body axis. This reduces the thrust component due to the angle of attack 
dependence that would appear in the energy rate equation. 
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Figure 1. Horizontal plane intercept geometry 
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Figure 2. Flight envelope for the F-8 aircraft 
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profiles for cases 1-3 



Figure 7. Commanded and actual altitude 
profiles for cases 1-3 
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Figure 10. Comparison of full boundary layer optimization with 
sub-optimal proportional lift controller 
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